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DESCRIPTION  OF  MECHANICAL  BEHAVIOR  OF  INELASTIC  SOLIDS 


E.  To 


Onat 


SUMMARY 

The  paper  is  concerned  with  the  representation  of  mechanical 
behavior  of  non-linear  solids  with  memory.  The  nature  and  limi¬ 
tations  of  integral  representation  of  Frechet  are  discussed. 
Definitions  of  state  and  state  variables,  which  are  based  on  obser¬ 
vable  histories,  are  introduced.  These  notions  give  rise  to  dif¬ 
ferential  equation  representation  of  behavior  discussed  in  the  latter 
part  of  the  paper. 


Most  of  the  results  in  this  paper  were  obtained  in  the  course 
of  research  sponsored  by  the  Department  of  Defense,  Advanced  Re¬ 
search  Projects  Agency,  under  supervision  of  Ballistic  Research 
Laboratories,  Aberdeen  Proving  Ground,  under  Contract  DA-19-020-AMC- 
0077(R),  ARPA  Order  No .71,  with  Brown  University. 
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Professor  of  Engineering,  Brown  University.  Now  at  Department  of 
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1.  INTRODUCTION 
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In  elastic  solids  the  stress  depends  on  the  deformation  existing  at  the  time 
of  measurement  of  the  stress.  For  inelastic  solids,  such  as  viscoelastic  or  plastic 
materials,  this  assumption  is  no  longer  valid.  The  current  stress  in  such  a  solid 
may  depend  in  some  manner  and  degree  upon  the  deformations  which  existed  in  it  prior 
to  the  time  of  measurement.  From  the  point  of  view  of  classical  continuum  mechanics 
the  description  of  mechanical  behavior  of  inelastic  solids  would  therefore  involve 
functional  relationships  between  the  stress  and  the  deformation  gradient  histories. 
Green  and  Rivlin  [1]  and  later  Green,  Rivlin  and  Spencer  [2]  have  shown  that  the 
dependence  of  the  stress  on  the  history  of  deformation  gradient  cannot  be  arbitrary. 
It  is  subject  to  two  major  restrictions.  The  first  of  these  arises  from  the  fact 
that  a  simultaneous  rotation  of  the  deformed  body  and  reference  system  must  leave 
the  stress  components  unaltered.  The  second  restriction  comes  from  any  symmetry 
which  the  material  may  possess  in  its  reference  state.  These  authors  have  also  con¬ 
structed  integral  representations  of  the  functional  dependence  of  the  stress  on  the 
history  of  deformation  gradient.  A  recent  review  of  modern  developments  in  the  con¬ 
tinuum  mechanics  of  inelastic  solids  has  been  given  by  Rivlin  [3]. 

The  present  paper  is  concerned  with  a  critical  review  of  the  subject  from  the 
point  of  view  of  an  experimenter.  The  above  mentioned  theoretical  work  provides  the 
experimenter  with  the  form  or  structure  of  the  relationship  which  exists  between 
the  histories  of  stress  and  strain.  The  task  of  filling  this  structure  with  phys¬ 
ical-  information  is  left  to  the  experimenter.  It  is  also  up  to  the  experimenter  to 
discover  whether  a  given  structure  constitutes  a  convenient  tool  for  the  description 
, of  mechanical  behavior  of  a  given  solid.  It  is  hoped  that  the  present  survey  can 
be' of  help  to  the  experimenter  in  the  performance  of  this  double  task. 

For  reasons  of  simplicity,  most  of  the  paper  is  restricted  to  the  study  of 
a  one-dimensional  situation.  Two  brief  introductory  sections  are  followed  by  a 
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section  on  the  strain  response  to  a  simple  family  of  stress  histories.  This  section 
serves  to  illustrate  a  fundamental  step  in  the  study  of  functionals. 

It  is  shown  in  Sections  4  and  5  that  for  certain  solids  Frechet’s  integral 
representation  of  mechanical  behavior  would  involve  a  prohibitively  large  number  of 
multiple  integrals  for  the  stress  or  strain  histories  of  interest.  A  more  con¬ 
venient  representation  of  mechanical  behavior  for  such  solids  could  be  developed 
by  generalizing  the  notion  of  piece-wise  smooth  functions  to  functionals.  An 
attempt  in  this  direction  was  recently  made  by  Onat  and  Wang  [4]. 

Sections  7  and  8  are  devoted  to  a  discussion  of  the  differential  equation 
representation  of  mechanical  behavior.  This  mode  of  representation,  which  is 
intimately  connected  with  the  notions  of  state  and  state  variables,  has  not  as 
yet  received  the  full  attention  of  the  workers  in  the  field,  especially  in  the 
area  of  finite  deformations. 

Section  7  contains  a  brief  discussion  of  the  notions  of  state  and  state 
variables  and  illustrates  the  circumstances  which  allow  one  to  represent  a  given 
history  dependence  by  a  system  of  differential  equations.  It  is  noted  in  this 
section  that  the  differential  equations  representation  constitutes  a  natural  tool 
for  the  description  of  the  piece-wise  smooth  behavior  encountered  in  plastic  or 
nearly  plastic  solids. 

Ordinary  plastic  solids,  such  as  structural  metals  at  room : temperature  and 
at  very  slow  strain  rates,  are  time  independent;  the  current  strain  in  such  ma¬ 
terials  depends  on  the  path  of  stress  history  in  the  stress  space,  but  it  is  in¬ 
dependent  of  the  speed  with  which  the  stress  point  moves  on  this  path.  In  spite 
of  this  simplification  it  is  not  convenient  to  express  the  path  dependence  of 
the  strain  by  means  of  an  integral  representation  because  of  lack  of  analyticity 
associated  with  the  question  of  loading  and  unloading.  In  structural  metals  sub¬ 
jected  to  high  strain  rates  or  temperatures  the  current  strain  depends  on  both 
the  stress  path  and  the  speed  of  the  stress  point  on  this  path  and  the  above 
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mentioned  difficulty  associated  with  loading  and  unloading  continues  to  exist.  The 
representation  of  mechanical  behavior  is  more  difficult  in  this  case  than  in 
ordinary  plastic  solids;  however  the  differential  equation  representation  for  plas¬ 
tic  solids  which  is  discussed  briefly  at  the  end  of  the  paper  can  easily  be  modi¬ 
fied  to  account  for  the  rate  effects. 

The  last  section  of  the  paper  is  concerned  with  the  questions  of  repre¬ 
sentation  in  the  case  of  arbitrary  deformations. 

2  o  PHENOMENOLOGICAL  STUDY  OF  MECHANICAL  BEHAVIOR 


We  wish  to  review  briefly  the  main  steps  involved  in  the  phenomenological 
study  of  mechanical  behavior  of  a  solid.  For  reasons  of  simplicity  we  first 
confine  ourselves  to  the  simple  case  of  tensile  test  which  enables  one,  in  prin¬ 
ciple,  to  study  homogeneous  deformations  of  an  initially  isotropic  solid  under 
the  time  dependent  uniaxial  tension.  We  shall  be  interested  in  small  isothermal 
deformations. 

The  study  begins  with  thq  preparation  of  a  number  of  specimens  which  at 
1  !*« 

the  start  of  testing  are  identical  in  every  respect, 

A  typical  experiment  involves,  if  a  stress  controlling  machine  is  at  our 

disposal,  application  of  time  dependent  stress  c(t)  on  the  time  interval  [0,T] 

and  observation  of  the  resulting  time  dependent  strain  e(x)  on  the  same  interval, 
:  T  T 


The  result  of  such  an  experiment  is  an  input-output  pair  (o(i),  e(t))  or  for 

:  •  o  o 

short  (a,e)  on  the  observation  interval  [0,T].  Different  inputs,  say,  o^\ 


(n) 


may  be  applied  to  different  identical . specimens  to  obtain  corresponding 


outputs  e 


(1) 


(n) 


The  result  of  such  a  multiple  experiment  is  a  set  of 


input -output  pairs  (q^?-,e^V)  . 


(a 


(n)  (n) 


). 


Reproducibility  of  tests  is  a  necessary  condition  for  the  specimens  to  be  iden¬ 
tical.  Sufficient  conditions  for  identity  are  more  difficult,  if  not  impossible, 
to  elucidate. 
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The  experimenter  would  consider  that  a  sufficient  insight  is  gained  into  the 
mechanical  behavior  of  the  solid  by  a  multiple  experiment,  if  for  the  stress  his¬ 
tories  (loading  programs)  a(f)  of  interest  he  can  guess  with  reasonable  accuracy 
what  the  corresponding  strain  responses  e(t)  would  be  without  performing  further 
tests.  When  this  stage  is  reached  the  experimenter  would  like  to  translate  the 

i 

knowledge  on  mechanical  behavior  thus  gained  into  the  mathematical  language  for 
the  purposes  of  transmitting  it  to  the  designer  or  the  analyst. 

The  task  of  the  experimenter  is  then  to  discover  and  to  represent  the 

relationship  that  exists,  for  a  given  inelastic  solid,  between  stress  and  strain 
A 

histories. 

Problems  similar  to  the  one  discussed  above  arise  in  the  study  of  com¬ 
munication  and  control  systems  under  the  name  of  Characterization  and  Representa¬ 
tion  Problems  or  simply  the  Identification  Problem  [5].  It  is  hoped  that  the 
workers  in  mechanics  would  find  it  profitable  to  follow  the  considerable  amount 
of  progress  which  is  being  made  in  the  study  of  the  identification  problem. 

Before  closing  this  section  we  must  point  out  that  the  study  of  the  mech¬ 
anical  behavior  is  often  much  more  difficult  than  the  previous  remarks  would  in¬ 
dicate.  For  instance  in  tests  involving  short  times,  such  as  wave  propagation 
and  impact  experiments,  the  measurement  of  stresses  is  nearly  impossible,  so  that 
the  experiment  provides  one  with  wave  shapes  and  speeds  measured  at  various 
locations  instead  of  a  set  of  input-output  pairs  (cr,e).  The  task  of  determining 
mechanical  behavior  of  the  solid  from  the  latter  evidence  is  an  involved  one  but 
it  is  not  unrelated  to  the  class  of  problems  discussed  in  this  paper. 


* 

A  more  thorough  study  of  mechanical  behavior  would  involve  histories  of 
temperature,  heat  flow,  etc.  Inclusion  of  these  would  complicate  the  picture 
but  would  introduce  no  new  conceptual  difficulties. 
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3c  MATHEMATICAL  PRELIMINARIES 

We  now  introduce  several  basic  concepts  which  are  indispensable  to  the  study 
of  input-output  pairs.  We  do  this  not  in  the  abstract  but  by  referring  to  the 
physical  situation  discussed  above. 

We  start  by  assuming  that  at  time  x  =  0  there  is  available  a  set  of  identical 
specimens  composed  of  the  given  solid.  The  specimens  carry  no  stress  at  this  time 
and  the  strains  for  t  £  0  are  measured  with  respect  to  the  configuration  at  t  =  0 
so  that 

a(x)  =  0,  e(x)  =  0  when  x  =  0.  (1) 

We  restrict  the  attention  to  the  loading  programs  a(x)  which  are  continuous 
on  the  observation  interval  [0,T]  and  vanish,  as  (1)  indicates,  at  x  =  0.  We  shall 
assume  that  the  strain  responses  e(x)  corresponding  to  continuous  strain  histories 
o(r)  are  also  continuous. 

In  each  test  the  solid  assigns  to  a  given  a(x)  a  response  e(x)  on  the  inter¬ 
val  [0,T].  An  electrical  engineer  would,  therefore,  consider  the  solid  or  the 
specimen  as  a  black  box.  A  mathematician  would,  on  the  other  hand,  regard  the  solid 
as  an  operator  and  would  use  the  following  notation, 

e  =  F[c]  (2) 

to  express  the  relationship  between  continuous  input -output  pairs  defined  on  [0,T]. 

One  would  remain  outsi.de  the  realm  of  mathematics  or  physics  if  one  did  not 
assign  some  properties  of  regularity  to  the  operator  F.  Here  we  shall  assume  that 
the  operator  F  for  the  solid  of  interest,  or  for  short,  the  solid  F,  is  continuous. 
For  the  purposes  of  the  present  paper  it  may  suffice  to  adopt  the  following  de¬ 
finition  of  continuity: 


We  further  assume  that  the  solid  is  in  the  zero  state  at  x  =  0;  this  means  that 
o(x)  =  0  for  x  £  0  implies  e(x)  =  0  for  x  *  0'.:  1 


r-  6  - 

It  is  said  that  the  operator  F  is  continuous ,  if  for  each  e  >  0  there  exists 
a  6  >  0  such  that  the  inequality 

jE(1)(x)  -  e<2)(,T)|  <  e,  on  [0,T]  (3) 

holds  for  the  strain  responses  produced  by  any  two  stress  histories  which . satisfy 
the  condition 

|o(1)(t)  -  a(2)(x)|  <  6,  on  [0,T].  (4) 

In  more  physical  terms  continuity  of  F  implies  that  the  strain  responses 
obtained  in  tests  conducted  with  slightly  different  loading  programs  differ  from 

j'c 

each  other  only  slightly. 

it  may  be  safe  to  assert  that  most  solids  of  interest  exhibit  this  property. 
However,  as  is  well  known,  the  Voigt  solid  constitutes  an  important  exception.  In 
this  solid  slightly  different  strain  histories  may  produce  different  stress  his¬ 
tories. 

Another  important  property  of  F  follows  from  the  observation  that  the 
happenings  in  the  future  cannot  affect  the  present  and  therefore  the  strain  at 
time  t  depends  only  on  the  stresses  on  the  interval  [0,t].  This  observation  im¬ 
plies  that  the  operator  F  must  have  the  property  of  causality  expressed  by  the 
symbolism' 

t 

e(t )  =  F[a( t )] .  •  (5) 

o 

The  above  equation  represents,  for  a  given  t,  a  functional  of  the  stress  histories 
defined  on  [0,t]. 

Another  property  which  we  may  assign  to  F  is  related  to  the  question  of 
"aging" .  The  state  of  the  solid  at  t  =  0  may  be  such  that ,  under  the  absence  of 

it 

Here  the  word"diff erent" is  used  in  the  sense  of  inequalities  in  (3)  and  (4). 
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stress  and  strain  for  t  z  0,  mechanical  behavior  of  the  solid  may  not  be  affected 
by  the  passage  of  time:  that  is  to  say,  if  e(x)  and  a(f)  is  an  input-output  pair, 
then  the  following  histories  would  also  constitute  an  input -output  pair: 

fo  ( 0  o  <  t  <  x 

e^f)  ,  cjl(t)  =<  (6) 

/  e(t-x)  I  a(t-x)  x  <  t 

for  any  non-negative  x.  Solids  which  have  the  above  property  are  called  time- 
invariant  or  non-aging  solids.  It  can  be  shown  that  for  time -invariant  solids 
the  operator  F  may  be  replaced  once  and  for  all  by  a  single  functional  defined 
on  the  interval  [0,T], 

4.  STRAIN  RESPONSE  TO  SIMPLE  LOADING  PROGRAMS 


In  this  section  we  wish  to  illustrate  several  concepts  introduced  in  the 
previous  sections.  The  present  section  will  also  prepare  the  ground  for  a  dis¬ 
cussion  of  integral  representations  of  mechanical  behavior. 

/We  consider  the  set  of  loading  programs  (Fig.  1)  defined  by 


c(x)  = 


T"tl 

'°i  +  ‘vV  t^t 


2  1 


0  £  T  £  t . 


\  *  t  <:  t2 


(7) 


where  t^  and  t^  are  given  fixed  times  and  and  are  arbitrary  constants 
satisfying  the  inequalities 


0  -  and  a2  S  M.  (8) 

We  wish  now  to  focus  our  attention  on  the  dependence  of  the  strain  e  at  time  t^ 
on  the  stress  histories  of  the  type  (7),  Since  each  loading  program  belonging 
to  the  above  set  is  fully  characterized  by  the  magnitudes  and  a 2,  the  continu¬ 
ous  functional  F  in  (5)  reduces,  for  argument  functions  defined  in  (7),  to  a 
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continuous  function  of  and  a 

_t 

e(t)  =  F[o£t)']'  =  f(a1,o2)„  (9) 

The  above  passage  from  continuous  functionals  to  continuous  functions  constitutes 
a  key  step  in  the  study  of  functionals. 

We  must  emphasize  again  that  in  (7)  t^  and  t2  are  kept  constant  but  and 
are  varied  so  that  the  dependence  of  f  on  t^  and  t^  does  not  enter  into  the  present 
cons iderat ions . 

A  multiple  experiment,  composed  of  n  inputs  of  the  type  (7)  and  the  corres¬ 
ponding  outputs,  determines  the  value  of  f  at  n  points  in  (8).  If  these  points  are 

chosen  suitably  then  the  experimenter  would  obtain  a  very  good  idea  as  to  what  sort 

A 

of  function  f  is;  he  would  then  be  able  to  construct  an  interpolating  function  f 
which  would  constitute  an  approximate  representation  of  f  even  over  the  points  which 
are  not  covered  by  the  multiple  experiment.  We  note  that  the  choice  of  n-points 
may  depend  on  the  nature  of  f  (for  instance,  more  experimental  points  would  be  needed 
over  the  subregions  of  (8)  where  f  is  "steep").  The  mathematical  representation  of 

A 

the  interpolating  functions  f  would  also  depend  on  the  nature  of  f,  as  we  shall 
presently  see. 

We  wish  now  to  review  the  form  of  f  for  some  basic  solids.  For  linear 

elastic  solids  we  have,  of  course, 

a2 

e(t)  =  f(clSa2)  =  -g-  .  (10) 

On  the  other  hand,  for  linear  viscoelastic  solids 

f  =  a^  +  a2a2’  (H) 

where  a^  and  a2  are  constants.  We  note  that  in  this  case  f,  the  present  value  of 
strain,  depends  on  the  stress  applied  at  previous  instances  through  the  presence 
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of  c^  in  (11).  For  most  polymers  (11)  constitutes  an  excellent  approximation  for  f 
provided  that  c^  and  and  hence  M,  is  sufficiently  small.  For  higher  values  of  M, 
f  becomes  non-linear.  In  some  cases  the  non-linearity  in  f  can  be  adequately  repre¬ 
sented  by  adding  higher  order  times  to  (11).  For  instance  an  experimental  study  of 
polypropylene  fibers  by  Ward  and  Onat  [63  indicated  that  up  to  the  extensions  of  2% 


f 


2 

a^,  +  a„a0  +  £  a..,  a. a.  a.  . 

11  2  2  ij,k  Vk  M  * 


(12) 


The  adequacy  of  the  above  representation  may  indicate  that  the  partial  derivatives 
of  f  up  to  the  fourth  order  are  continuous  in  the  domain  of  interest. 

We  next  consider  the  rigid-plastic  material  defined  by  the  stress-strain 
diagram  shown  in  Fig.  2  where  Y  is  the  yield  stress  and  is  the  tangent  modulus 
which  for  reasons  of  simplicity  is  assumed  to  be  constant.  The  function  f(a1,o2), 
which  describes  the  mechanical  behavior  of  this  material  for  stress  histories  (7) 
is  therefore  given  by 


(13) 


Note  that  although  f  is  continuous,  its  first  derivatives  exhibit  discontinuities 
on  the  interfaces  of  the  domains  indicated  in  (13). 

The  following  remarks  concerning  f  defined  in  (13)  will  be  useful  in  the 
discussion  of  integral  representations.  Since  f  is  continuous  it  can  be  approxi¬ 
mated  to  any  degree  of  accuracy  by  a  polynomial  P  according  to  the  theorem 

of  Weierstrass  (see  [7],  p.  481).  However  the  degree  of  the  Weierstrass  polynomial 
which  achieves  an  acceptable  approximation  over  a  closed  region  of  interest  may  be 
very  high  due  to  the  presence  of  discontinuities  in  3f/3o^  and  3f/3o2  in  (13). 

To  elaborate  on  this  remark  we  may  consider  the  simple  case  of  Y  =  0  and  =  1 


where  f  becomes: 
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f  =  max(a  ,a2)  (14) 

We  now  consider  the  polynomial  approximation  characterized  by  the  inequality  and 

e  ft 

the  domain  given  below  : 

|max(CT1,o2)  -  P(ol5o2)|  <  e  on  0  <  cr1  ,a2  £  1  (15) 

It  can  be  shown,  using  a  result  of  Bernstein  [8]  concerning  the  related  function 
y  =  |x|  that  the  degree  N  of  the  polynomial  P(a^,a2)  in  (15)  must  at  least  be 
N  £  »  This  result  indicates  once  more  that  the  Weierstrass  polynomial  is  not 

a  natural  tool  for  the  representation  of  the  piece-wise  linear  function  defined 
in  ( 13 )  „ 

5.  INTEGRAL  REPRESENTATION  OF  MECHANICAL  BEHAVIOR 


We  begin  with  linear  solids.  A  solid  is  said  to  be  linear  if  the  operator 
F  representing  the  mechanical  behavior  of  the  solid  has  the  properties  of  homo¬ 
geneity  and  additivity: 


(a)  F[Ao]  =  AF[o]  where  A  is  a  constant. 


(b)  F[a(1)+  c(2)]  =  F[c(1)]  +  F[a(2)], 


The  last  property  implies  that  the  strain  response  of  the  material  to  the  sum 
of  two  stress  histories  is  equal  to  the  sum  of  the  strain  responses  to  the  indi¬ 
vidual  stress  histories  (the  principle  of  superposition). 

Riesz's  well-known  work  on  linear  functionals  [9]  when  applied  to  the  pre¬ 
sent  case  provides,  for  differentiable  stress  histories,  the  following  integral 
representation: 


It  will  be  assumed  throughout  the  paper  that  all  relevant  quantities  have  been 
properly  nori-dimensionalized. 


e(t)  = 


J(t,x) 


(16) 


4  o 

where  the  function  J(t,x)  is  independent  of  the  argument  function  and  characterizes 
the  linear  solid. 

When  the  solid  is  time -invariant  then  J  becomes  a  function  of  the  single 
argument  (t-x).  The  types  of  experiments  needed  for  the  establishment  of  linearity 
of  a  given  solid  and  for  the  measurement  of  the  associated  creep  function  J(t) 
or  the  relaxation  modulus  G(t)  are  well-known  and  will  not  be  discussed  here. 

Many  solids  of  interest,  such  as  metals  subjected  to  high  stresses,  tempera¬ 
tures  or  stress  rates,  exhibit  inelastic  behavior  but  do  not  possess  the  properties 
(a)  and  (b)  of  linearity. 

Representation  of  mechanical  behavior  of  such  nonlinear  solids  is  a  more 
difficult  problem.  We  have,  of  course,  the  following  result  based  on  Frechet's 
work  on  nonlinear  functionals  [11]  which  is  at  our  disposal  for  the  purposes  of 
representation : 

AA 

A  continuous  and  causal  operator  F  acting  on  the  set  of  equicontinuous 
functions  a(x)  on  [0,T]  can  be  represented  to  any  desired  degree  of  accuracy  by 
a  sum  of  multiple  integrals: 


t 

F  CcJ(x)] 
o 


■t 

K1(t,x1)a(x1)dx1  +  . . . 
o 

ft  ft 


(17) 


+ 


Km(t’V 


x  )a(x.  ) . . .a(x  )dx  _ dx 

m  l  ml  m 


o  o 

where  the  kernels  are  determined  by  the  operator  independently  of  the  argu¬ 
ment  function  o(x)  and  they  are  continuous  functions  of  x^.  In  the  case  where  F 
is  time-invariant,  K.  become  continuous  functions  of  t-x.. 

’  l  l 


For  a  review  of  this  subject  see  [10], 

AA 

For  a  definition  of  equicontipuity  or  compactness  see  [12]. 
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If  we  are  to  use  (17)  to  represent  mechanical  behavior  of  a  given  solid  then 
we  must  conduct  a  multiple  experiment  aimed  at  the  determination  of  (a)  the  number 
of  terms  to  be  retained  in  (17)  for  a  desired  degree  of  accuracy  and  (b)  the  kernels 
defining  these  terms . 

A  multiple  experiment  based  on  the  loading  programs  (7)  discussed  in  the 
previous  section  would  shed  light  on  both  items.  We  observe  first  that  for  stress 
histories  (7)  the  sum  (17)  would  become  a  polynomial  of  degree  m  in  and 
If  the  experiments  suggest  that  can  be  represented  adequately  by  a 

polynomial,  say,  of  third  degree  then  the  indication  would  be  that  the  first  three 
integrals  should  be  retained  in  (17).  One  could  then  conduct  further  tests  in¬ 
volving  more  complicated  strain  histories,  say,  histories  composed  of  three  or 
four  linear  segments,  to  see  whether  this  initial  guess  remains  a  good  one. 

Such  tests  would  also  help  to  determine  the  structure  of  the  relevant  kernels  K^. 

Studies  which  follow  this  approach  have  been  conducted  by  Ward  and  Onat  [6], 
Lifshitz  [13]  and  Onaran  and  Findley  [14]  for  non-linear  polymers. 

Indications  are  that  for  polymers  studied  in  the  above-mentioned  work,  non¬ 
linearity  in  behavior,  which  occurs  when  stresses  become  relatively  large,  can  be 
adequately  represented  by  the  first  few  terms  of  the  sum  (17). 

For  metals  a  different  situation  appears  to  exist.  We  have  shown  in  the 
previous  section  that  for  a  rigid  plastic  material  the  representation  of 
by  a  polynomial  would  require  a  prohibitively  larger  number  of  terms.  This  would 
mean  therefore  that  a  large  number  of  terms  would  be  needed  in  the  Frechet  sum 
for  the  same  material.  At  first  sight  this  would  not  cause  any  concern  since  other 
means  of  representation  are  known  for  such  solids.  However,  recent  work  by  Wang 

and  Onat  [15]  on  creep  of  an  alumihum  alloy  at  elevated  temperatures,  showed  that 

'  * 

for  this  material  the  representation  of  the  strain  response  £(0^,02)  to  loading 
programs  composed  of  two  steps  of  magnitudes  0^  and  0 ^  would  also  involve  a  poly¬ 
nomial  of  high  degree  in  0^  and  0 This  is  partly  due  to  the  fact  that  the 
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function  z(o^,o^)  behaves  differently  over  the  domains  characterized  by  the  in¬ 
equalities  and  a  >  a^. 

This  observation  indicates  that  aluminum  retains  the  piece-wise  smooth 
behavior  associated  with  plasticity  in  spite  of  the  occurrence  of  viscous  effects 
at  high  temperature. 

The  above  remarks  show  that  there  exist  solids  of  interest  for  which  the 
Frechet  description  of  mechanical  behavior  would  involve  a  large  number  of  terms 
for  stress  histories  of  interest .  Such  solids  may  be  said  to  be  strongly  non¬ 
linear.  Since  a  constitutive  law  of  the  type  (17)  containing  a  large  number  of 
terms  is  likely  to  be  cumbersome  in  applications,  it  may  be  appropriate  to  look  for 
more  convenient  representations  of  mechanical  behavior  for  strongly  nonlinear 
solids. 

For  this  purpose  one  could  try  to  generalize  the  notion  of  piece-wise 
smooth  functions  to  operators  and  functionals.  Representation  of  a  piece-wise 
smooth  functional  would  involve  a  sum  of  the  type  (17)  containing  a  few  multiple 
integrals  whenever  the  argument  functions  lie  in  a  given  sub-domain  of  the  space 
■of  argument  functions.  When  one  goes  to  an  adjacent  sub-domain  then  a  different 
but  equally  simple  integral  representation  would  take  over.  An  attempt  in  this 
direction  has  been  made  by  Onat  and  Wang  in  [4],  where  they  discuss  order  sensi¬ 
tive,  piece-wise  smooth  (or  linear)  functionals. 

Representation  of  mechanical  behavior  could,  in  certain  cases,  be  based 
on  a  system  of  differential  equations.  As  will  be  seen  in  sections  6  and  7,  this 
mode  of  representation  is  intimately  connected  with  the  notions  of  state  and  of 
state  variables.  We  will  see  in  section  8  that  representation  by  differential 
equations  is  particularly  suitable  for  strongly  nonlinear  solids. 
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6=  REPRESENTATION  OF  MECHANICAL  BEHAVIOR  BY  DIFFERENTIAL  EQUATIONS  IN  LINEAR 
VISCOELASTICITY 

In  this  introductory  section  and  in  the  sections  which  will  follow,  it  will 
be  more  convenient  to  regard  strain  histories  as  inputs  and  stress  histories  as 
outputs. 

We  consider  a  time-invariant  linear  solid  defined  by  the  constitutive  law: 

“t 

c(t)  =  G(t-x)e(x )dx .  (18) 

o 

We  assume  that  the  relaxation  modulus  G(t)  is  composed  of  a  sum  of  exponentials: 

n  -t/T. 

G(t )  =  I  G.e  1,  (19) 

i=l  1 

where  G^  and  T^  are  positive  constants. 

We  introduce  the  following  n  quantities 

*  -(t-x )/T. 

q.  =  e  e°(x)dx,  (20) 

o 

and  observe  by  taking  the  time  derivative  of  (20)  that 

q^t) 

q!(t)  +  — -  =  e°(t),  (21) 

i 

On  the  other  hand  a  combination  of  (18),  (19)  and  (20)  provides  the  equation 

n 

o(t)  =  I  .G.q.(t).  (22) 

i=l  11 

We  now  observe  that  (21)  and  (22),  together  with  the  initial  conditions 

q^  =  0  when  t  =  0  ( 23 ) 

which  follow  from  (20)  can  be  regarded  as  the  differential  equation  representation 
of  the  mechanical  behavior  of  the  solid  (18). 

It  is  important  to  note  from  (21)  and  (22)  known  at  a  given 

■  da:  QRhtjftD, 
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time  t,  then  a(t)  can  be  determined  for  x  >  t  from  the  knowledge  of  e'(x)  or 

equivalently  from  the  knowledge  of  e(x)  -  e(t)  for  t  >  t.  Thus  if  the  strains 

are  measured  for  x  i  t  with  respect  to  the  configuration  at  time  t,  the  behavior 

of  the  solid  on  [t,T]  will  depend  only  on  the  strains  on  [t,T]  and  on  the  at 

time  to  This  remark  suggests  that  q^  may  be  called  state  variables.  The  existence 

of  a  finite  number  of  state  variables  implies  that  the  mechanical  behavior  of 

the  solid  does  not  depend  on  all  the  details  of  the  strain  history.  As  can  be 

seen  from  (20),  those  "components"  of  the  strain  rate  histories  e‘(x)  which  are 

x/T. 

"orthogonal"  to  the  functions  e  do  not  contribute  to  q^  and  therefore  do  not 

ft 

affect  the  future  behavior  of  the  material. 

7.  THE  CONCEPT  OF  STATE.  REPRESENTATION  BY  DIFFERENTIAL  EQUATIONS. 

Consider  a  set  of  specimens  which  have  been  subjected  to  some  strain 
histories  on  the  interval  [0,t].  We  apply  to  these  specimens  further  strains 
on  Ct.Tj.  We  denote  the  strain  on  [t,T]  based  on  the  configuration  at  time  t 

A 

by  e  (x).  In  the  present  case  of  small  strains  we  have,  of  course, 

e(x)=e(x)-e(t),  x  £  t  (24) 

where  e(t)  and  e(x)  are  the  strains  referred  to  the  zero  state. 

We  shall  say  that  a  set  of  specimens  are  in  the  same  state  at  time  t  if 
any  pair  of  identical  strain  histories  e  (x)  on  [t,T]  produce  in  any  two  of  these 
specimens  identical  stress  histories  on  [t,T]. 

The  test  pair  shown  in  Fig,  3  suggests  that  the  strain  histories  and 

(2) 

e  on  [0,t]  may  have  produced  the  same  state  at  time  t.  In  order  to  make  sure 

that  this  is  actually  the  case  one  must  conduct  further  test  pairs  using  the  same 

(1)  (2)  ft 

strain  histories  e  and  e  on  [0,t],  but  different  e  (x)  on  [t,T]„ 

ft 

See  [16]  for  a  further  discussion-  of  this  topic. 
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A  strain  history  on  [0,t]  produces  a  given  state  at  time  to  It  will  be  of 
interest  to  know  the  number  of  states  at  time  t  created  by  all  strain  histories  on 
[0,t]»  This  number  will,  in  general,  be  infinite  and  some  of  our  future  remarks 
will  deal  with  comparison  of  infinities. 

For  an  elastic  solid  all  strain  histories  on  [0,t]  which  have  the  same 
terminal  value  e(t)  produce  the  same  state  at  time  t.  Thus  all  possible  states 
of  an  elastic  solid  at  time  t  can  be  represented,  in  the  present  uniaxial  case, 
by  points  on  the  real  axis,  each  point  having  the  coordinate  e(t).  We  may 
therefore  say  that  all  possible  states  of  an  elastic  solid  fill  the  real  axis 
(or  a  segment  of  it,  if  the  magnitude  of  strains  of  interest  is  bounded). 

Next  we  consider  a  Maxwell  solid  in  its  familiar  model  representation 
(Fig.  4 ) o  It  can  easily  be  seen  from  the  definition  of  the  state  that  two  strain 
histories,  which  at  time  t  give  rise  to  the  same  elongation  in  the  spring  of  the 
model,  produce  the  same  state.  Therefore  all  possible  states  of  the  Maxwell  model 
also  fill  the  real  axis. 

On  the  other  hand,  for  the  linear  viscoelastic  solid  considered  in  section  6, 
a  state  is  characterized,  as  we  have  seen  before,  by  the  n  quantities  so  that 
the  state  space,  for  this  material,  is  n-dimensional.  The  extreme  case  occurs  when 
every  distinct  strain-history  on  [0,t]  produces  a  different  state.  In  this  case 
the  state  space  coincides  with  the  space  of  input-functions. 

We  now  restrict  the  discussion  to  materials  which  possess  an  n-dimensional 
state  space.  The  coordinates  of  the  state  point  will  be  denoted  by  q^.  During 
a  test  the  state  point  will  move  in  the  state  space  starting  from  the  point 
corresponding  to  the  state  at  t  =  0. 

We  wish  now  to  prove  that  the  stress  at  time  t  must  be  a  function  of  q^(t) 

and  1 1 

o(t)  =  f(q^(t) ,t) . 


(25) 
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The  proof  follows  from  the  definition  of  state  and  from  the  fact  that  the  operator 
F  defining  the  mechanical  behavior  of  the  solid  assigns  continuous  a(x)  to  continu¬ 
ous  e(x)  on  [0,T]o 

Consider  two  strain  histories  on  [0,T].  Assume  that  at  time  t 

these  histories  give  rise  to  the  same  state.  The  following  continuous  strain 
history 


e(x) 


|  e^(x)  on  [0,t] 
l^e(1)(t)  +  e5£(2)(i)  on  [t,T] 


(26) 


will  produce  on  [t,T],  by  the  definition  of  the  state,  the  same  stress  history 
as  that  produced  by  But  the  stresses  caused  at  time  t  by  e^Vand  e 

defined  in  (26)  must  be  the  same  by  the  above  mentioned  property  of  the  operator 
F,  Therefore  the  stresses  at  time  t  produced  by  and  e^2^  must  be  the  same, 

which  proves  the  assertion. 

It  can  also  be  shown  that  if  the  solid  is  time -invariant  then  o  in  (25) 
will  depend  only  on  q.(t). 

We  now  consider  the  dependence  of  the  state  variables  q^  on  the  history  of 
strain.  It  follows  from  the  definition  of  state  that 


A 

qi  =  ( t )  a t ; e  ]  on  [t,T] 


(27) 


The  above  symbol  means  that  q^(x)  on  [t,T]  is  related  by  the  operator  F^  to  the 

A 

history  of  strain  e  (t)  based  on  the  configuration  at  time  t,  the  operator  itself 
being  a  function  of  q^(t)  and  t.  We  may  assume  that  the  operator  F^  is  continuous 
and  causal. 

Now  we  observe  that  q^  becomes  a  given  function  of  x  on  [t,T]  for  a  given 
& 

history  e  (x)„  It  is  of  interest  to  study  the  derivative 


_d_ 

dx 


qi<T)X: 


q:(t) 


(28) 
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It  can  be  shown  that  if  the  operator  possesses  a  Frdchet  derivative  and  if 
q^(x )  for  e  (x)  =  0  on  [t,T]  has  a  derivative  at  x  =  t,  then 

q°(t)  =  g.(q.(t),t)  +  h.(q. (t),t)e°(t),  (29) 

3-  1  J  •*-  J 

where  and  Ik  are  functions  of  the  indicated  arguments. 

Equations  (25)  and  (29)  together : with  the  initial  conditions  q^(0)  =  q? 
constitute  a  differential  equation  representation  of  the  mechanical  behavior 
of  the  solid. 

That  e"  appears  linearly  in  (29)  is  due  to  the  existence  of  the  Frechet 
derivative  of  F^  which  may  be  too  strong  a  requirement  to  place  upon  F^  for  some 
solids.  It  is  likely  that  for  solids  such  as  metals  where  the  behavior  depends 
strongly  on  whether  there  has  been  loading  or  unloading  in  the  recent  past  (29) 
will  have  the  more  general  form 

q° (t )  =  f .(q. ,e*  ,t)  (30) 

±  ±  J 

and  f^  will  not  necessarily  be  continuous  in  its  arguments. 

All  that  has  been  said  would  remain  a  rather  formal  mathematical  structure 
if  one  cannot  give  broad  rules  concerning  the  types  of  experiments  needed  for  the 
determination  of  the  number  of  state  variables,  and  of  the  functions  governing 
their  growth.  It  would  seem  that  there  exists  no  systematic  study  of  such 

is 

questions.  It  should  be  remembered  however  that  one  rarely  encounters  an 
entirely  unknown  material.  With  each  class  of  materials  there  is  available  a 
certain  amount  of  knowledge  of  phenomenological  or  microscopic  nature  which  pro¬ 
vides  starting  points  in  such  studies.  A  brief  study  of  plasticity  in  the  next 
section  will  illustrate  this  point. 

* 

Rules  based  on  simple  histories  of  the  type  (7)  should  not  be  difficult  to 
construct . 
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8 „  ARBITRARY  DEFORMATIONS 


Green  and  Rivlin  [1]  have  shown  that  the  functional  dependence  of  stress 

& 

on  the  deformation  gradient  history  must  be  of  the  form 

t  m 

g(t)  =  e(t)  F  [e£r)]  D1 (t )  (31) 

where  g  is  the  stress  tensor,  D(t)  is  the  tensor  of  deformation  gradients  at 
T 

time  t,  0  (t)  is  its  transpose  and  F  is  a  tensor  valued  functional  of  the  strain 
history  on  [0,t].  Here  the  tensor  g(t)  is  an  appropriate  measure  of  finite  strain. 
The  above  form  is  such  that  a  simultaneous  rotation  of  the  deformed  body  and 
reference  frame  leaves  the  stress  component  unaltered. 

In  the  case  where  the  time-invariant  solid  is  continuous  and  the  components 
e..(f)  of  the  strain  history  belong  to  a  set  of  equicontinuous  functions  a  straight - 
forward  generalization  of  (17)  provides  the  following  integral  representation 
for  F : 


=  K.  .  (t-T)e,  .  (r)dx  + 
i]  o  kl  kl 


+  Kijklmn(t-tl’t-T2)Ekl(Tl)i:n,n(T2)dTldT2  +  — 


If  the  solid  possesses  symmetry  in  its  reference  state  at  t  =  0,  then  the  tensor 
valued  kernels  in  the  above  expression  must  exhibit  certain  properties  of  invari¬ 
ance.  These  have  been  studied  exhaustively  in  recent  years  (see,  for  instance, 
Rivlin  [3]). 

In  the  particular  case  of  isotropic  solids  the  above  expression  takes  the 
following  remarkable  form  [1]: 


We  assume  that  at  t  =  0  elements  of  the  solid  are  identical  and  are  identically 
oriented  with  respect  to  a  fixed  coordinate  frame. 
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F  S  K  I  +  E 
°  B=1 


Kg(t-T 1 ,  t-Tg) 


(33) 


Cg(t1)  ..o  eC t g )  +  s(fg)  ...  e(x1)]dT1  ...  dx^ 


where  I  is  the  unit  tensor  and 

ft  ft 

<j>CX(t-x1, 

0 


s  are  polynomials  in  the  invariants  1^  given  by 
....  t-xB)  tr  gCxp  ...  g(xg) 


(34) 


dx1dx.2 


dx. 


these  polynomials  being  dependent  on  the  arguments  shown  in  ( 33 ) . 

If  the  strains  are  small  enough  so  that  terms  of  degree,  say,  higher  than 
the  third  may  be  neglected  in  comparison  with  those  of  the  third  degree,  the 
expression  (33)  can  be  simplified  considerably  [17].  It  can  easily  be  seen  that 
the  simplified  expression  contains  twelve  unknown  functions.  Lockett  [17]  has 
shown  how  these  functions  could,  in  principle,  be  determined  experimentally  by 
means  of  stress-relaxation  experiments.  Experimental  studies  based  on  this  ap¬ 
proach  have  been  reported  by  Lifshitz  [13]  and  Onaran  and  Findley  [14], 

The  success  of  this  approach  depends  on  the  range  of  validity  of  the 
simplified  expression.  As  remarked  before.,  for  certain  polymers  the  range  of 
validity  may  be  adequately  large,  but  for  metals  it  may  be  so  small  that  a 
different  approach  to  the  problem  of  representation  may  be  needed. 

In  the  discussion  of  the  differential  equation  representation  of 
mechanical  behavior  we  restrict  ourselves  to  small  deformations.  For  small 
deformations  (small  shape  changes,  small  rotations)  (31)  reduces  to 

t 

o(t)  =  F[e(x)]  (35) 

o  * 

where  g  now  is  the  strain  tensor  employed  in  the  classical  theory  of  elasticity. 
If  one  considers  a  time -invariant  material  equations  (25)  and  (30)  take 
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the  form 


°ij(t)  =  fij(<Ja(t>)' 

(36) 

qa(t)  =  fa(qe(t)’ 

where  and  are*  respectively,  the  components  of  the  stress  and  strain 
tensors,  in  a  fixed  rectangular  frame. 

An  important  question  which  must  be  faced  at  this  stage  is  concerned  with 
the  possible  tensorial  character  of  the  state  variables.  The  state  variables  can 
be  regarded  as  average  quantities  describing  the  microscopic  structure  of  the 
solid  and  therefore  it  would  be  natural  to  expect  that  they  will  be  tensor  valued. 

It  can  also  be  argued  that  it  would  be  more  convenient  to  work  with  tensor 
valued  state  variables.  For  instance,  if  the  state  variables  are  tensor  valued 
then  one  could  apply  recent  results  [18]  on  the  form  invariance  to  (36)  to  obtain 
a  representation  appropriate  to  isotropic  solids. 

ieit 

We  now  consider  briefly  the  representation  of  elastic-plastic  behavior. 

It  is  convenient  in  this  case  to  regard  the  stress  histories  as  inputs.  Constitu¬ 
tive  equations  for  an  elastic-plastic  solid  have  the  following  structure,  which 
is  similar  to  (36): 


ei:(t)  *  cijkiaki(t)  +  fij(q«<t)>’ 


fakl(0ii(t)'Vt))0kl(t) 


The  first  term  in  (37)  represents  elastic  components  of  the  strain,  the  second 
represents  the  plastic  components.  The  latter  depend  on  the  history  of  stress 


Generalization  of  the  definitions  of  state  and  state  variables  to  the  present 
case  of  arbitrary  small  deformations  is  self-evident  and  will  not  be  given  here. 

For  a  detailed  discussion  of  the  stress-strain  relations  in  Plasticity,  see  [19]. 


22 


and  hence  are  functions  of  the  state  variables  q  . 

In  (38)  the  case  (a)  occurs  when  the  current  stress  point  is  on  the  yield 
surface  and  is  moving  out  of  it  at  the  instant  of  interest: 

S'F 

f(a..(t),  q  (t))  =  0  and  -r - a!.  >0.  (a) 

ID  a  30^  i] 

The  function  f  which  defines  the  yield  surface  depends  on  the  history  through 

the  state  variables  q  . 

a 

The  case  (b)  occurs  when  (i)  the  stress  point  is  within  the  yield  surface 
or  (ii)  the  stress  point  is  on  the  yield  surface,  but  it  is  not  moving  out: 


f  <  0,  or  f  =  0  and  -r - a’..  <0.  (b) 

’  So.,  id 

ID 

It  may  be  of  interest  to  comment  here  on  the  number  n  of  state  variables. 

In  a  similar  treatment  of  Plasticity,  Kroner  [20]  uses,  as  state  variables,  the 
components  of  a  symmetric  second  order  tensor  which  is  a  measure  of  the  disloca¬ 
tion  loop  density.  Green  and  Naghdi's  treatment  of  Plasticity  [21]  involves,  in 
the  present  context,  seven  state  variables,  six  of  them  being  the  components  of 
the  plastic  strain  tensor. 

In  his  theory  of  strain-hardening,  Prager  [22]  assumes  that  the  yield  sur¬ 
face  moves  in  the  stress  space  with  no  change  in  shape  and  orientation.  Therefore 
Prager* s  theory  involves  six  state  variables. 

Experiments  are,  of  course,  the  ultimate  source  of  information  on  the 
number  of  state  variables.  A  multiple  experiment  based  on  different  stress  his¬ 


tories  will  produce  a  family  of  yield  surfaces.  The  number  of  parameters  needed 
for  a  (approximate)  description  of  this  family  can  be  taken,  as  a  first  guess, 


as  the  number  of  state  variables. 
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